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Abstract 

Homodyne tomography provides a way for measuring generic field-operators. Here 
we analyze the determination of the most relevant quantities: intensity, field, ampli- 
tude and phase. We show that tomographic measurements are affected by additional 
noise in comparison with the direct detection of each observable by itself. The case 
of of coherent states has been analyzed in details and earlier estimations of tomo- 
graphic precision are critically discussed. 



1 Introduction 



One of the most exciting developments in the recent history of quantum optics 
is represented by the so-called Homodyne Tomography, namely the homodyne 
detection of a nearly single-mode radiation field while scanning the phase of 
the local oscillator [1-4]. ^From a tomographic data record, in fact, the density 
matrix elements can be recovered, thus leading to a complete characterization 
of the quantum state of the field. This is true also when not fully efficient pho- 
todetectors are involved in the measurement, provided that quantum efficiency 
is larger than the threshold value rj = 1/2. 

In homodyne tomography a general matrix element is obtained as an expec- 
tation value over homodyne outcomes at different phases. In formula 



7T +00 



l^lv 9 ) = J ' — J dxp v (x; 0) (i/)\K v (x - x^ip) 



where p v (x] cf>) is the probability density of the homodyne outcome x at phase 
<p for quantum efficiency rj and the integral kernel is given by 
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K v (x) = -Re y d/c/c exp I -k 2 + ikxj . (2) 



o 



While the kernel in Eq. (2) is not even a tempered distribution, its matrix 
elements can be bounded functions depending on the value of 77. This is the 
case of the number representation of the density matrix, for which the "pattern 
function" 



f^fafl^inlK^x-xJln + d) . (3) 

can be expressed as a finite linear combination of parabolic cylinder functions 
[3]. 



As it comes from the experimental average 



& = J — j dxp v (x; 0) fi v) n+d (x, 0) , (4) 
-00 

the tomographic determination g n ^ m for the matrix element Q n ,m = (n\g\m) 
is meaningful only when its confidence interval is specified. This is defined, 
according to the central limit theorem, as the rms value rescaled by the number 
iV of data. As Q n ^ m is a complex number, we need to specify two errors, one 
for the real part and one for the imaginary part respectively. For the real part 
one has 



1 /2 
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where 
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Km = / V / tePiM) [Re/2 +d (*,0)] 2 > ( 6 ) 

-00 

and likewise for the imaginary part. 

Quantum tomography opened a fascinating perspective: in fact, there is the 
possibility of device-independent measurements of any field-operator, includ- 
ing the case of generalized observables that do not correspond to selfadjoint 
operators as, for example, the complex field amplitude and the phase. The first 
application in this direction has been presented in Ref. [5] where the number 
and the phase distributions of a low excited coherent state have been recovered 
from the original tomographic data record. No error estimation was reported 
in Ref. [5], whereas an earlier analysis of the precision of such determinations 
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has been reported in Ref. [6] on the basis of numerical simulations. The idea 
behind these papers is simple. Any field operator A, in fact, is described by its 
matrix elements A n , rn = (n\A\m) in the number representation. Then, upon a 
suitable truncation of the Hilbert space dimension, at the maximum photon 
number H, the expectation value of A is given by the linear combination 



H 

(A) = Qn,m A n ^ m , (7) 

n,m=0 

whereas the corresponding confidence interval is evaluated by error propaga- 
tion calculus 



H 
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AA ~ J2 \ A Qn,m lAvnl ■ (8) 
n,m=0 

The whole procedure relies on two assumption, namely 

n , m <C 1 n, m > H , (9) 

and 



lim A£ n , m = 0, (10) 



n,m— > oo 



which needs a more careful analysis. The condition in Eq. (9) is certainly ful- 
filled for some value of H, whose determination, however, requires an a priori 
knowledge of the state under examination. On the other hand, it has been 
shown in Ref. [7] that in a tomographic measurement involving N experimen- 



tal data the errors Ree„ iJn and lme n:tn saturate to the value y2/iV for 77 = 1, 
whereas they diverge exponentially for r\ < 1. Therefore, the condition (10) 
cannot be fulfilled in a real experiment and we conclude that determinations 
of Ref. [5] are not meaningful, as they are affected by diverging errors and are 
based on a priori knowledge of the state. For the same reason the analysis of 
Ref. [6] is not correct, and the added noise has been largely overestimated. 



2 Homodyning field operators 



In this paper we analyze the tomographic determination of field-quantities 
from a different perspective. By homodyning an observable A we mean the 
average 
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of the state-independent kernel function 7?.L4](:r;0) [8], which allows for the 
determination of the expectation value (A) without the detour into density 
matrix elements. For a Hilbert-Schmidt operator A Eq.(ll) follows directly 
from a generalization of Eq.(l) with 7?.L4](ie; 0) = Tr ^AK(x — £</,)}, whereas 
alternative approaches to derive explicit expressions of the kernel have been 
suggested [10,8], that here we briefly recall. Starting from the identity involving 
trilinear products of Hermite polynomials (valid for k + m + n = 2s even [9] ) 



Jdxe^H k (x) H m (x) H n (x) = {s _ ms r m) T {s n : n)] , (12) 

— oo 

Richter proved the following nontrivial formula for the expectation value of 
the normally ordered field operators [10] 



N ' J tj: J ./2 n + m ( n+m \ 

-co v \ n J 

which corresponds to the kernel 

n[a^ n a m ](x; 0) = e i(m-n)<t> J^±M^) (13) 

For nonunit quantum efficiency the homodyne photocurrent is rescaled by 
7] whereas the normally ordered expectation (a) n a m ) gets an extra factor 
772 ("+m)_ Therefore, one has 



7e>tn a ™ ](x . 0) = e i (ro -n)» f^f^J , (14) 

V(2^) n+m ( t ) 

where the kernel TZ V [O] (x; 0) is defined as in Eq. (11), but now with the 
experimental probability distribution p v (x;<j)) for nonunit quantum efficiency 
77. From Eq. (14) by linearity one can obtain the kernel lZ v [f](x; 0) for any 
operator function / that admits a normal ordered expansion 



/ = /(a,at)= £ /S^V"- (15) 

ram=0 



4 



One obtains 



n+m,. 



E 



nm=0 

H s (^jx)i s d s 



j s \(2rj) s / 2 dv' 



v=0 



(16) 



where 



nm=0 



m 



(17) 



Continuing from Eq. (16) one obtains 



K v [f}(x] 4>) = exp 
and finally 



2rj dv 2 dv 



(18) 



v=0 



+oo 

TZ v [f](x; 0) = J -^= e -l w2 F[f}(w + 2ix; 0) . (19) 

— oo 

In summary, the operator / possesses a tomographic kernel TZ v [f](x; 0) if the 
function F[f\(v\ 0) in Eq. (17) grows slower than exp(— r/v 2 /2) for v — > oo. In 
addition, as we can assume that p v (x] 0) goes to zero faster than exponentially 
at x — > oo, the average in Eq. (11) is meaningful for the integral in Eq. 
(19) growing at most exponentially for x — > oo. In the next section we will 
consider the tomographic determination of four relevant field quantities: the 
field intensity, the real field or quadrature, the complex field, and the phase, 
for all of which the above conditions are satisfied. 



3 Added noise in tomographic measurements 



As already mentioned in the previous section the tomographic measurement of 
the quantity A is defined as the average of the kernel w v = 1Z V [A] (x, 0) over 
the homodyne data. A convenient measure for the precision of the measure- 
ment is given by the confidence interval Aw, which, being w v a real quantity, 

is given by Aw v = (w 2 — w^ 2 \ , where 
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w 2 = n 2 {A}(x,<j ) ) = J -f Jdx Pri (x,<j>)1ll[A](x,<j>). (20) 

-oo 

When the quantity A can also be directly measured by a specific setup it makes 
sense to compare tomographic precision Aw with the corresponding fluctua- 
tions \J (AA 2 ). Notice that, when we deal with r) < 1 the noise \J (AA 2 ) V is 
larger that the quantum fluctuations due to smearing effect of nonunit quan- 
tum efficiency. As we will see, the tomographic measurement is always more 
noisy than the corresponding direct measurement for any observable, and any 
quantum efficiency rj. However, this is not surprising, in view of the larger 
amount of information retrieved in the tomographic measurement compared 
to the direct measurement of a single quantity. 

In Table 1 we report the tomographic quantities w v the field-observables exam- 
ined. Before going into details of each observable, we mention a useful formula 
for evaluating confidence intervals. These are obtained by averaging quantities 
like 



7^[a tn a m ](x,0) = e 2 ^ (m 



H 2 n+m (y/Ejx) 



(2^) ( n + m ) 

By means of the following identity for the Hermite polynomials [11] 



2 • 



(21) 



m + n 
m 



*w = ^ £ «,;?■[-*)■ ■ (22) 



we arrive at 



|2 |2 m+n (Oh\\n k 
t?L H 'VV ^m+n k\ A (n + m-k)\ »U A > W» V V 

which expresses the generic square kernel 7^[a' i ' ri a m ](x, 0) in terms of "diago- 
nal" kernels TZ v [a^ h a k ](x, 0) only. 



5. 1 Field- Intensity 



Photodetection is the direct measurement of the field-intensity. For a single- 
mode of the radiation field it corresponds to the number operator n = a^a. 
For nonunit quantum efficiency rj at the photodetectors, only a fraction of the 
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incoming photons is revealed, and the probability of detecting m photons is 
given by the Bernoulli convolution 



oo 

Pv( m ) = 2 P™ 



>n 



n=m 



r) m (l - r])" , (24) 



m 



p nn being the actual photon number distribution of the mode under examina- 
tion. One considers the reduced photocurrent 



I = -a) a , (25) 
V 

which is the quantity that traces the photon number, namely it has the same 
mean value 



A) = - Yl m p( m ) = n ' ( 26 ) 

^ m=0 

where we introduced the shorthand notation n = (a' a). On the other hand, 
one has 



— 1 00 /l \ 

{AP) V = — J2 m 2 p{m) = (An 2 ) + n - - 1 . (27) 

V m=0 K 1 ! J 

In Eq. (27) (An 2 ) denotes the intrinsic photon number variance . The term 
n{r]~ l — 1) represents the noise introduced by inefficient detection. The tomo- 
graphic kernel that traces the photon number is given by the phase-inde- 
pendent function w v = 2x 2 — (2rf)~ 1 . With the help of Eq. (23) we can easily 
evaluate its variance, namely 



2_3\ J_ 



Aw, = (An 2 ) + -(n 2 ) +n ( ---)+_. (28) 



The difference between Aw v and (AI 2 ) V defines the noise N[h] added by 
tomographic method in the determination of the field intensity 



N[n] = - 



(n 2 ) + n(--l) +\ 



(29) 



The noise N[n] added by the tomographic measurement is always a positive 
quantity and largely depends on state under examinations. For coherent states 
we consider the noise-ratio 
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Aw v * r i / i s 

dn v = a — = < 2 + - r/n H 



1/2 



(30) 



which is minimum for n = i] 1 . 



3.2 Real Field 



For a single mode light-beam the electric field is proportional to a field quadra- 
ture x = 1/2 (g^ + a), which is just traced by homodyne detection at fixed 
zero-phase with respect to the local oscillator. The tomographic kernel, that 
traces the mean value Tr{^£} is given by w v = TZ v [x](x,<f)) = 2xcos0, inde- 
pendently on r], whereas the square kernel w 2 = 7Z 2 [x](x, 0) = Ax 2 cos 2 <fi can 
be rewritten as 



w 2 = j [}Z[a 2 ](x, 4>) + n[aV](x, 0)] + 1Z[Ja](x, <f>) + i- . 



The confidence interval is thus given by 



(31) 



Aw r , 



(a^) + (a 2 ) 



1 / a + 

+ n H ( 

2 V \ 2 



(A^> + in + ^,(32) 



(Ax 2 ) being the intrinsic quadrature fluctuations. For coherent states Eq. (32) 
reduces to 



Ay; 



n + 



(33) 



The tomographic noise in Eq. (32) has to be compared with the rms variance 
of homodyne detection for nonunit quantum efficiency This is given by 



(Ax 2 ), = (Ax 2 ) + 



1—7] 

4:7] ' 



(34) 



For coherent states Eq. (34) becomes (Ax 2 ) r/ = 1/(47/) The added noise results 



N[x] 



(35) 



whereas the noise-ratio for coherent states is given by 
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and increases with the scaled intensity rjn. 



3.3 Field amplitude 



The detection of the complex field amplitude of a single-mode light-beam 
is represented by the generalized measurement of the annihilation operator 
a. The tomographic kernel for a is given by the complex function w v = 
lZ[a\(x,(f)) = 2xexp(i(f>). To evaluate the precision of the measurement one 
has to consider the noise of a complex random variable. Generally there are 
two noises 



Aw r 



— \w ri \ 2 ± | Ait;^| 



(37) 



corresponding to the eigenvalues of the covariance matrix. Using Eq. (23) one 
has 



w v = K v [a\(x,(f)) = e 



i2<f> 



- + 2n„[a^a](x,<f)) 
V 



and 



K| 2 = \K v [a](x,(f)\' 



7/ 



which lead to 



D i2<t> 



+ n ri {a 2 ](x,<i)) . (38) 



1 + 2r]Tl v [a t a](x,(j)) 



(39) 



Aw r 



I + 2 n-|(a)| 2 ± 



(a 2 ) - (a)'' 



(40) 



because e m< ^ = S n Q for all states. The optimal measurement of the complex 
field a, corresponding to the joint measurement of any pair of conjugated 
quadratues x^ and X0 +7r / 2 can be accomplished in a number of different ways: 
by heterodyne detection [12], eight-port homodyne detection [13-15], or by 
six-port homodyne detection [16,17]. In such devices each experimental event 
a = x + iy in the complex plane consists of a simultaneous detection of the 
two commuting photocurrents x and y, which in turn trace the pair of field- 
quadratures. The probability distribution is represented by the generalized 
Wigner function W s (a,a) with ordering parameter s related to the quantum 
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efficiency as s = 1 — 2i] 1 . The precision of such measurement is defined 
likewise Eq. (37) as follows 



(Aa 2 ), 



1 



lap — lal 2 ± 



a 2 — a 2 



(41) 



where 



a — J d 2 a aW s (a, a) = (a) 
c 

a 2 = J d 2 a a 2 W s (a, a) = (a 2 ) 
c 

/I 
d 2 a aa* W s (a, a) = (a^a) H — 
1 1 



(42) 



^From Eqs. (41) and (42) we have 



(Aa 2 ), 



n+--\(a)\ 2 ± 



[a 2 ) - (a) 2 



(43) 



The noise added by quantum tomography is thus simply given by 



N[a] = -n , 

which is independent on quantum efficiency. 
For a coherent state we have 



(44) 



1 

n H — 
rj 



and (Aa 2 ), = — , 



and the noise ratio is given by 



(45) 



5a r 



\ (Aa 2 ), 



{1 + i]n} 



1/2 



(46) 



3.4 Phase 



The canonical description of the quantum optical phase is given by the prob- 
ability operator measure [18,19] 
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d/i((f)) — — ^ exp{i(m — n)0}|n) (m\ , (47) 

2 71 " n ,m=0 



which defines a phase operator [20] through the relation 



1 

}= f diitf) = -i £ (-) n_m \n){m\ . (48) 

In principle, a comparison between homodyne tomography and direct deter- 
mination of the phase would require from one side the average of the kernel 
corresponding to the operator 0, and from the other side the direct experi- 
mental sample of the operator 0. However, such a comparison would be purely 
academic, as there is no feasible setup achieving the optimal measurement (47). 
For this reason, here we consider the heterodyne measurement of the phase, 
and compare it with the phase of the tomographic kernel for the correspond- 
ing field operator a, i.e. w v = aig(2xe^). Notice that the phase w v is not just 
the given local oscillator phase, because x has varying sign. Hence averaging 
w v is not just the trivial average over the scanning phase 0. The probability 
distribution of such kernel variable can be easily obtained by the following 
identity 

TT OO TT ^ OO 

J — J dxp v (x,(f)) = 1 = J Jdxp v (x,w v ) , (49) 

~oo —IT 

which implies 



p v {w v ) = - dxp v (x,w v ) . (50) 

7T J 


The precision in the tomographic phase measurement is given by the rms 
variance Aw^ of the probability (50). In the case of a coherent state \p) = | 
(zero mean phase) Eq. (50) becomes 



Pv( w v) = ^ j 1 + Erf 

which approaches a "boxed" distribution in [— 7r/2, tx/2] for large intensity. We 
compare the tomographic phase measurement with its heterodyne detection, 
namely the phase of the detected complex field a. The outcome probability 
distribution is the marginal distribution of the generalized Wigner function 
W s (a, a) (s = 1 — 2r/~ 1 ) integrated over the radius 



V2\P\ 



cosw. 



n 



(51) 
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p„(0) = JpdpW s (pe l tpe-^), 
o 



(52) 



whereas the precision in the phase measurement is given by its rms variance 
Af v . We are not able to give a closed formula for the added noise N[(f>] = 

2 2 

Aw v — A<f> . However, for high excited coherent states \p) = (zero mean 

phase) one has Ay v = tt 2 /12 and A(/> = (2rjn)~ 1 . The asymptotic noise-ratio 
is thus given by 



5<p v = 




-V w > 1 



(53) 



A comparison for low excited coherent states can be performed numerically. 
The noise ratio 5^ (expressed in dB) is shown in Fig. 1 for some values of 
the quantum efficiency r). It is apparent that the tomographic determination 
of the phase is more noisy than the heterodyne one also in this low-intensity 
regime. 



4 Summary and Remarks 



Homodyne tomography provides a complete characterization of the state of 
the field. By averaging suitable kernel functions it is possible to recover the 
mean value of essentially any desidered field-operator. In this paper we ana- 
lyzed the determination of the most relevant observables: intensity, real and 
complex field, phase. We have shown that these determinations are affected 
by noise, which is always larger than the corresponding one from the direct 
detection of the considered observables. In Table 2 a synthesis of our results 
is reported. 

We have considered the ratio between the tomographic and the direct-measu- 
rement noises. This is an increasing function of mean photon number n, how- 
ever scaled by the quantum efficiency rj. Therefore homodyne tomography 
turns out to be a very robust detection scheme for low quantum efficiency. In 
Fig. 2 the noises ratio (in dB) for all the considered quantities are plotted for 
unit quantum efficiency versus n: this plot has to compared with Fig. 6 of Ref. 
[6], where the tomographic errors were largely overestimated. 
In conclusion, homodyne tomography adds larger noise for highly excited 
states, however, it is not too noisy in the quantum regime of low n. It is then a 
matter of convenience to choose between a direct measurement and homodyne 
tomography, as the former is the most precise measurement of the desidered 
quantity, whereas the latter represents the best compromise between the con- 
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flicting requirements of a precise and complete measurement of the state of 
radiation. 
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Table 1 



Tomographic vs. direct quantities for the variables of interest in this paper. 
VARIABLE TOMOGRAPHIC QUANTITY DIRECT QUANTITY 



Intensity 


w v 


= 2x 2 - -i- 


I = 


- a^a 


Real Field 


w v 


= 2x cos 4> 


x = 


= 5(« + « f ) 


Complex Amplitude 


w v 


= 2x exp{ 


a = 


= x + iy 


Phase 


w v 


= arg(xe^) 




= arg(a) 



Table 2 

Added noise in tomographic determinations and noise ratio for coherent states. For 
the phase the results are valid in the asymptotic regime n S> 1 

VARIABLE ADDED NOISE NOISE RATIO 



Intensity 


N[n] 


l 

~~ 2 


y> + ft (a _ i) + ^" 


5n v 




Real Field 


N[x] 


1 

— 2 






5x v 


= {2(l+>,ft)} 1/2 


Complex Amplitude 


N[a] 


= \ n 




Sa v 


= {1 + r/n} 1 ^ 2 


Phase 


N[cf>] 


IT 

~ 12 


i_ 

2-qn 
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Fig. 1. Ratio between tomographic and heterodyne noises in the measurement of 
the phase for low excited coherent states, The noise ration is reported versus the 
mean photon number n for some values of the quantum efnency. From bottom to 
top we have 7] = 0.2, 0.4, 0.6, 0.8, 1.0. 
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Fig. 2. The noises ratio (in dB) for all the quantities considered in this paper. They 
are plotted for unit quantum efficiency versus n: this plot has to compared with 
Fig. 6 of Ref. [6]. 
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